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ABSTRACT 

Background Imaging of Cosmic Extragalactic Polarization (BiCEP) is a bolometric polarimeter designed to 
measure the inflationary Z?-mode polarization of the cosmic microwave background (CMB) at degree angular 
scales. During three seasons of observing at the South Pole (2006 through 2008), BiCEP mapped ^ 2% of the 
sky chosen to be uniquely clean of polarized foreground emission. Here, we present initial results derived from 
a subset of the data acquired during the first two years. We present maps of temperature. Stokes Q and U, E and 
B modes, and associated angular power spectra. We demonstrate that the polarization data are self-consistent 
by performing a series of jackknife tests. We study potential systematic errors in detail and show that they are 
sub-dominant to the statistical errors. We measure the ii-mode angular power spectrum with high precision 
at 21 < i < 335, detecting for the first time the peak expected at £ ^ 140. The measured E-mode spectrum 
is consistent with expectations from a ACDM model, and the B-mode spectrum is consistent with zero. The 
tensor-to-scalar ratio derived from the B-mode spectrum is r = 0.02^11^, or r < 0.72 at 95% confidence, the first 
meaningful constraint on the inflationary gravitational wave background to come directly from CMB Z?-mode 
polarization. 

Subject headings: cosmic background radiation — cosmology: observations — gravitational waves — infla- 
tion — polaiization 



1. EVITRODUCTION 

One of the cornerstones in our current understanding of cos- 
mology is the theory of inflation. Inflation addresses several 
major shortcomings of the standard big bang model, resolving 
the flatness and horizon problems and explaining the origin of 
structure; however, the theory has yet to be unambiguously 
confirmed by observational evidence. 

Numerous experiments have demonstrated that the cosmic 
microwave background (CMB) is an extremely effective tool 
for studying the early universe. Precision measurements of 
the temperat ure anisotropies now span a wide range of an- 
gulai- scales (Jones et al."2006'i Reichar dt et"aD l2009'; 'Nolta' 
enil. 2009; Friedman et al. 2009, Sieve rs et al.|2 009; Brown 
et al.||2009| and have yielded tight constraints on a model of 
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the universe in which the energy content is dominated by a 
cosmological constant and cold dark matter (ACDM). 

The polarization anisotropies of the CMB provide even 
more insight into the history of the universe, potentially en- 
coding information from long before the moment of matter- 
radiation decoupling. The primary source of CMB polar- 
ization is Thomson scattering of the local quadrupole of the 
photon-baryon fluid sourced by density fluctuations. The re- 
sulting partial polarization has no handedness and is called the 
gradient or "/i-mode" by analogy to curl-free electric fields. 
The ACDM parameters that predict the temperature spectrum 
also predict the shape of the E-mode spectrum with almost 
no additional information (reionization enhances power at the 
largest an gular scales). E-mo de polarization was first detected 
by Dasi (Kovac et al.|20 02[) and has since been measured by 
many other expe riments (Heitch et al. 2005! 'Montroy et aL| 
,2006,|Sievers et al. 2007 , Wu et al. 2007 ; Bischoff et al. 2008 
The acoustic peaks in the polarization spectra have been mea- 
sured to high precision in bo th the TE spectrum ( [Piacentiiu] 
et al.|2005 Nolta et al.|2009|l and, more recently, direc tly in 
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the EE spectrum (Pr yke et al.|2009[ [Brown et al |20091 l, pro 
viding further support for our basic understanding of CMB 
physics. 

Inflation predicts the existence of a stochastic gravitational 
wave background, created during the initial accelerated ex- 
pansion of the universe, which imparts a unique imprint on 
the CMB at the surface of last scattering (Polnarev 1985). 
In addition to producing E-mode polarization, gravitational 
waves also induce a curl or "B-mode" in t he polarization 
anisot ropies ( Seljak & Zaldarriaga 1997 ; Ka mionkowski et aL| 
|1997| l. The inflationary B-mode signal is expected to peak at 
degree angular scales (multipole moment £ ^ 100) with an 
amplitude determined by the energy scale of inflation. Be- 
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cause density fluctuations at the surface of last scattering cre- 
ate only E-mode polarization, a detection of the B-mode sig- 
nal would b e strong evidence that inflation occurred (see, e.g., 
IDodelsone tal. 2009). 

The inflationary Z?-mode amplitude is parameterized by the 
tensor-to-scalar ratio r, and the most restrictive published 
upper limit, r < 0.22 (95% confidence), comes from mea- 
surements of large-scale temperature anisotropies in combi- 
nation wit h baryon acoustic os cillation and Type la super- 
nova data ( |Komatsu et al.|[2009( ). However, the constraints 
from temperature anisotropies are ultimately limited by cos- 
mic variance, and lowering the r limit further requires direct 
polarization measurements. Currently, limits from polariza- 
tion are still far worse than those from temperature — for ex- 
ample, assuming ACDM parameters that are fixed at Wmap 
best-fit values, the Wmap B-mode spectrum constrains r < 6 
at 95% confidence. The results reported in this paper provide 
upper limits on the B-mode signal that are an order of magni- 
tude more stringent than those set by Wmap. 

Background Imaging of Cosmic Extragalactic Polarization 
(BiCEP) is a microwave polarimeter that has been designed 
specifically to probe the B-mode of CMB polarization at de- 
gree angular scales. The instrument observed from the South 
Pole between 2006 January and 2008 December A detailed 
description of the BiCEP instru ment characterization p roce- 
dures is given in a separate paper, Takahashi et al. (2010 1, that 
accompanies this text. In this paper, we report initial CMB 
polarization results from the 2006 and 2007 observing sea- 
sons. 

2. THE BICEP INSTRUMENT 

A complete description of the BiCEP instrument is available 
Yoon et aL] ([2006 ), and only a brief summary is given here. 



in 



The BiCEP receiver consists of a two-lens refracting telescope 
coupled to a focal plane of 4 9 orthogonal pairs of polarization- 
sensitive bolometers (PSBs; pones et al.j2 003 ). The PSB pairs 
are divided between 25 that observe at 100 GHz and 24 at 
150 GHz (two of the 150 GHz PSB pairs were reconfigured 
for 220 GHz operation in late 2006 and were subsequently 
not used for CMB analysis). The angular resolution at 100 
and 150 GHz is 0.93° and 0.60°, respectively, and the instan- 
taneous field of view is 18°. The entire focal plane and op- 
tics assembly is housed in an upward-looking cryostat with 
toroidal liquid nitrogen and liquid helium tanks. The clean 
optical path and azimuthal symmetry minimize instrumental 
polarization systematics. 

The receiver is supported in an azimuth-elevation mount 
with a third degree of rotational freedom about the boresight. 
The mount is located on the top floor of the Dark Sector Lab- 
oratory (89.99° S, 44.65° W) at the Amundsen-Scott South 
Pole station, a site with excellent atmospheric transparency 
and stability at millimeter wavelengths as well as outstanding 
infrastructure. The telescope penetrates through the roof and 
is sealed to the building with a flexible environmental enclo- 
sure, leaving most of the instrument accessible in a warm lab 
setting. 

The 24-hr visibility of the target field from the South Pole 
enables uninterrupted observation and deep integration. Bl- 
CEP's primary CMB field lies within the "Southern Hole," a 
region of low dust emission outlined in Figure [1] in a right 
ascension and declination range of approximately \a\ < 60° 
and -70° < 6 < -45°. The telescope observation cycle is 48 
sidereal hours in length and is divided into four 9-hr CMB ob- 
servations, 6 hr of Galactic observations, and 6 hr of cryogen 
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Fig. 1. — BiCEP's CMB and Galactic fie lds are outlined on the 1 50-GHz 
FDS Model 8 prediction of dust emission ^Finkbeiner et al.|1999^ , plotted 
here in equatorial coordinates. 

servicing. The CMB field is covered twice over the same az- 
imuth range during each 48-hr cycle, but the elevation halves 
are mapped in opposite order between the two observations. 
The boresight angle in each cycle is held fixed at one of 
four angles {-45°, 0°, 180°, 135°} that provide good thermal- 
microphonic stability and redundant polarization angle cover- 
age. 

BiCEP maps the sky with azimuth-elevation raster scans. 
During each complete CMB observation (18 hr), the telescope 
boresight steps in elevation between 55° and 59.75° in 0.25° 
increments. At each step in elevation, the telescope performs 
a set of 50 back-and-forth azimuth scans over a total period of 
^50 minutes. The azimuth scan width is 64.4°, and the speed 
is held constant at 2.8°/s over ^70% of the scan duration, 
which modulates the sky signal and places it in a frequency 
band of approximately 0.1-1 Hz. The scans have a fixed az- 
imuth center that is updated at each elevation step to approx- 
imately track the field center This scan strategy was chosen 
instead of continuous tracking in order to allow removal of 
any azimuth-fixed or scan-synchronous contamination. 

Relative detector gains are measured regularly during ob- 
serving cycles with "elevation nods" performed at the begin- 
ning and end of each fixed-elevation scan set. During an 
elevation nod, the mount performs a rounded triangle wave 
motion in elevation with a 1.2° peak-to-peak amplitude, and 
the detector voltages vary in response to the changing line-of- 
sight air mass. The nod is performed over a 45-s period to 
reduce thermal disturbances on the focal plane, and thermal 
drifts are further suppressed by using mirror-image elevation 
nods at the beginning and end of each scan set (up-down-up 
and down-up-down). 

3. INSTRUMENT CHARACTERIZATION 

The timestream d(t) of a perfect linearly polarized detec- 
tor is related to the signal on the sky through the expression 
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d(t)= T{rh)+Q{rh)cos2'ip+U(rh)sm2'il), where the sky signal 
is described by the Stokes parameters T, Q, U. This expres- 
sion depends only on the detector's direction of observation Tb 
and polarization orientation angle i/;. (The time-dependence 
of and ip that arises from the scan strategy is suppressed 
for clarity.) Modifying this simple expression, the timestream 
response of a BiCEP PSB is described by 

c/(f) = ^(0*{n(0 + ^(0 J dvA,F(v) j drVir-Vb.v) 

[r(r , I/) + ^ (e(r , i.) cos 2 V + C/(r , u) sin 2^-)] } , ( 1 ) 
1 + e J 

where the extra terms are calibration quantities that describe 
the properties of the instrument. The cross-polar leakage, 
which is a PSB's level of response to orthogonally polarized 
light, is parameterized by e. The detector signal is convolved 
with the co-polar beam P(r, v), which has a spatial extent that 
depends on the coordinate r. The beam-convolved detector 
response is integrated over the pass band F{v), which is mul- 
tiplied by the effective antenna area A^. The gain factor g(t) 
converts voltage to temperature units, and n{t) is an additive 
noise term. Finally, the entire expression is convolved with 
the detector transfer function Kit). 

In order to faithfully reconstruct the temperature and po- 
larization signal on the sky, it is necessary to measure all the 
terms in Equation[T]that relate T, 2, U to the detector voltage. 
A complete description of these measurement procedures and 
results is giv en in the accompanyin g instrument characteri- 
zation paper (ITakahashi et al. 2010[ ). Here, we give only a 
brief summary of the quantities used directly in data analysis. 
The list includes detector transfer functions, absolute and rel- 
ative gains, main beam shapes, cross-polar leakages, detector 
orientatio n an gles, and pointing. The noise is discussed sepa- 
rately in §6.1 



3.1. Transfer Functions 

Deconvolving detector transfer functions is the first step 
in producing clean timestreams that are suitable for analysis. 
Relative gains are measured with elevation nods at 0.02 Hz, so 
the transfer functions must be characterized over a frequency 
range that spans at least 0.01-1 Hz in order to link the rela- 
tive gains to the entire science band. The transfer functions 
were measured with a microwave source (Gunn oscillator or 
broadband noise source) that was placed near the telescope 
window and square-wave modulated at 0.01 Hz. The time- 
domain responses to the transitions were Fourier transformed, 
divided by the transform of the modulation waveform, and av- 
eraged for each detector in order to obtain the deconvolution 
kernel. The transfer functions have a measurement precision 
of 0.5% rms across the signal band and have sufficiently high 
signal-to-noise ratio (S/N) to be deconvolved directly from 
the bolometer timestreams. (Although the transfer functions 
can be described by a simple model, we do not rely on those 
fits.) The relative gain uncertainty that results from transfer 
function deconvolution is < 0.3% over the range 0.01-1 Hz. 

3.2. Relative and Absolute Gains 

PSB relative gains are measured with elevation nods at the 
beginning and end of every set of constant-elevation azimuth 
scans. The nods inject an atmospheric signal into the bolome- 
ter timestreams, and the gains are obtained by fitting each 
timestream against the cosecant of the detector elevation. The 



resulting volts-per-air mass responsivity factors are normal- 
ized to the average of good detectors for each frequency band 
during the scan set. The common-mode rejection of CMB 
temperature anisotropics in gain-adjusted PSB pairs is mea- 
sured to be better than 98.9% at degree angular scales. 

We cross-correlate the CMB temperature fluctuations mea- 
sured by Wmap and BiCEP to obtain absolute gains, which 
relate CMB temperature to detector units. The Wmap maps 
are smoothed to BiCEP's resolution by applying the ratio of 
the beam window functions, Z?^'*-^^/Z?*^^^. (We do not ap- 
ply a correction for the pixelization of the Wmap maps since 
the effects are negligible.) The smoothed maps are then con- 
verted to simulated detector timestreams using the boresight 
pointing data. The timestreams are filtered and converted 
back into maps, thus creating a "BiCEP-observed" version of 
the Wmap data. The Bicep map and processed Wmap maps, 
which have compatible beam and filter functions, are cross- 
correlated in multipole space to obtain the absolute gain 



gb = 
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/ „WMAP-1 „* WMAP-2\ 



(2) 



where is a top hat binning operator, and are the spher- 
ical harmonic expansion coefficients. To avoid noise bias, the 
flfa, coefficients in the denominator are taken from two differ- 
ent Wmap maps; for this analysis, we have used the Q and 
V ban d maps from the five-year data release ( [Hinshaw et al.| 
|2009| l. The resulting gain calibration, is approximately 
flat over BiCEP's I range of 21-335, where the lower bound 
is set by the timestream filtering, and the upper bound is set by 
beam uncertainty. The absolute gain used for each of the Bi- 
CEP frequency bands is a single number taken from the aver- 
age of over six uniform bins spanning a multipole window 
of 56 < ^ < 265, and the absolute gain uncertainty is derived 
from the standard deviation of To assess the impact of 
errors in the beam window functions, we have calculated g\, 
using both the Q and V band Wmap data in the numerator 
of Equation |2] The average values are consistent within 
errors, and we take the larger of the two standard deviation 
values, 2%, as a conservative estimate of the absolute gain 
uncertainty in temperature units. 

3.3. Boresight and Detector Pointing 

The two components of detector pointing reconstruction are 
telescope boresight pointing and detector offsets relative to 
the boresight (focal plane coordinates). Raw boresight point- 
ing timestreams are obtained from encoders located on the 
three mount axes, and corrections to the raw data are applied 
from a model describing axis tilts and encoder offsets. The 
pointing model is established with star observations from an 
optical camera located on the upper surface of the cryostat. 

The focal plane coordinates of each detector are recon- 
structed from measurements of CMB temperature fluctua- 
tions. Temperature maps from each detector at each of the 
four boresight angles are cross-correlated with the full-season 
map, and the process is repeated with varying centroid adjust- 
ments until the maps are self-consistent. This method results 
in 0.03° rms centroid uncertainty. 

3.4. Polarization Orientation and Efficiency 

To recover polarization information from detector 
timestreams, it is necessary to know the polarization orien- 
tation angle ip and cross-polar leakage e of each PSB. Note 
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that e is a property of the detector itself and is independent 
of the cross-polar beam, which is a property of the optical 
chain upstream of the PSB. For BiCEP, e dominates over the 
cross-polar beam. Both tp and e are measured with several 
different devices (e.g., rotating wire grid, dielectric sheet 
calibrator) that send polarized light into BiCEP at many 
different angles with respect to the detectors. The phase of 
each PSB's sinusoidal response and the ratio of the minimum 
to maximum determine ip and e, respectively. The uncertainty 
in the measured orientation angles is ±0.7°. The median e for 
BiCEP PSBs is about 0.04, with a measurement uncertainty 
of ±0.01. 

3.5. Main Beam Shapes 

The BiCEP beams were mapped by raster scanning the tele- 
scope over a bright source at various fixed boresight angles. 
The beams are well described by a Gaussian model, with fit 
residuals typically about 1 % with respect to the beam ampli- 
tude. The beams at both frequencies are nearly circular, with 
ellipticities under 1.5%. We therefore approximate the beams 
as symmetric Gaussians with average full widths of 0.93° and 
0.60° at 100 and 150 GHz, respectively. The distribution of 
beam widths varies by ±3% across the focal plane, and each 
width is measured to a precision of ±0.5%. 

4. LOW-LEVEL ANALYSIS AND MAPMAKING 

The analysis presented here includes BiCEP data from the 
2006 and 2007 observing seasons. For these initial results, 
we restrict the data set to uninterrupted 9-hr CMB observa- 
tions taken during February-November Although there is 
no evidence for Sun contamination, we exclude data acquired 
during the Austral summer because of mediocre weather con- 
ditions and increased station activities. A significant fraction 
of 2006 was devoted to calibration measurements and inves- 
tigation of a different scan strategy, and the nominal data set 
in this season consists of 124 days. The amount of data in- 
creased to 245 days in the 2007 season. 

4.L Data Cuts 

The data set is further reduced by omitting 9-hr observation 
phases with extremely poor weather quality. For each phase, 
the standard deviation of relative gains from elevation nods 
is calculated for each channel, and the median over the 100 
and 150-GHz channels of those standard deviations yields two 
numbers per phase. An observation phase is cut if either of the 
median standard deviations is greater than 20% of the average 
relative gain. After this weather cut, the data set consists of 
1 17 and 226 days for the 2006 and 2007 seasons. 

Of the 49 optically active PSB pairs, several are excluded 
from each season due to anomalous behavior such as no re- 
sponse, excess noise, poorly behaved or time-dependent trans- 
fer functions, and exceptionally poor polarization efficiency. 
Six experimental pixels containing Faraday rotator modules 
(2006 season only) and the two 220-GHz pixels (2007 and 
2008 seasons) are also excluded from CMB analysis. A to- 
tal of 19 100-GHz and 14 150-GHz PSB pairs are used for 
2006 analysis, and 22 100-GHz and 15 150-GHz PSB pairs 
are used for 2007. 

4.2. Low-level Timestream Processing 

The raw output of BiCEP consists of voltage timestreams 
sampled at 50 Hz for 144 channels comprising 98 light 
bolometers, 12 dark bolometers, 16 thermistors, 10 resistors. 



and 8 intentionally open channels. The low-level timestream 
cleaning begins with concatenating and trimming the raw data 
files for each 9-hr observation phase. Complete transfer func- 
tions are deconvolved for the 98 light bolometers, and all 
timestreams are low pass filtered at 5 Hz and downsampled 
to 10 Hz. Relative detector gains are derived from eleva- 
tion nods, and horizon and celestial boresight coordinates are 
calculated using the pointing model. "Half-scans," or single 
sweeps in azimuth, are identified, and the turnarounds are ex- 
cluded from further analysis. The remaining central portion 
of the half-scan, which has nearly constant velocity, makes up 
~ 75% of the total scan duration. 

Each detector half-scan is subjected to three data quality 
checks. First, relative gains derived from elevation nods at 
the beginning and end of a set of scans are compared, and 
all of the half-scans in the scan set are excluded if the gains 
differ by more than 3%. Second, gain-adjusted PSB pairs are 
differenced over each half-scan, and the half-scan is rejected 
if the magnitude of the skew or kurtosis is abnormally high 
(PSB pair-difference timestreams are expected to be Gaussian 
white noise). Finally, half-scans that contain cosmic rays and 
other signal spikes are identified by points that are greater than 
7 times the standard deviation of the smoothed timestream. 
On average, the combined half-scan flagging criteria exclude 
about 3% of all half-scans over all light detectors. Because the 
flagged percentage is small, the problematic half-scans are not 
gap-filled and are simply omitted from analysis. 

4.3. Mapmaking 

After low-level cleaning, the bolometer timestreams are 
binned into temperature and polarization maps. We have de- 
veloped two data analysis pipelines for BiCEP that differ start- 
ing from the mapmaking stage. The code in each pipeline is 
completely independent of the other, and the only shared data 
products are the initial set of downsampled, cleaned detector 
timestreams, boresight pointing, and calibration data. Both 
pipelines have reproduced all the results reported in this pa- 
per, achieving similar statistical power and excellent agree- 
ment. The mapmaking algorithms used by the two pipelines 
are similar, although one (designated for this paper the "pri- 



mary pipeline") bins in the Healpix ( |G6rski et al.|2005] l 

1 for this pa- 



pixelization scheme, while the other (designated 
per the "alternate pipeline") produces maps in rectangular co- 
ordinates. In this section, we describe the primary pipeline's 
mapmaking p rocedure in detail. 

Following Jones et al. ( 2007 1, the simplified timestream 
output dij of a single PSB can be expressed as 

dij = gij [T(pj) + 7,(G(P>) cos 2^ij + U{T?j) sin l^Jij)] , (3) 

where g is the gain, T,Q,U are the beam-integrated Stokes 
parameters of the sky signal, 7 = (l-e)/(l+e)is the polar- 
ization efficiency factor, and ip is the PSB polarization ori- 
entation projected on the sky. The index / denotes the PSB 
channel number, j is the timestream sample number, and pj 
is the map pixel observed at time j. The goal of mapmaking 
is to recover T, Q, U from the bolometer timestreams. 

The mapmaking procedure for BiCEP begins with the for- 
mation of gain-adjusted sum and difference timestreams for 
each PSB pair: 



, 1 , 
dij = i^idiij / g2i.j±d2i+\.}/ g2i+\.j)- 



(4) 



To reduce atmospheric 1 // noise, a third order polynomial is 
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Fig. 2. — BiCEP T, Q, U. and coverage maps. The resolution is about 0.9° and 0.6° at 100 and 150 GHz, respectively, and no smoothing or apodizing has been 
applied to the maps. The noise per square degree in the central region of the Q and U maps is 0.81 /iK at 100 GHz and 0.64 fiK at 150 GHz. Note that the color 
scales of the temperature and polarization maps differ by a factor of 10. 
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subtracted from the sum and difference timestreams for each 
half-scan in azimuth. Azimuth-fixed and scan-synchronous 
contamination are removed by subtracting a template signal, 
which is formed by binning the polynomial-filtered detec- 
tor timestreams in azimuth over each set of fixed-elevation 
scans. There are slight differences in the scan-synchronous 
signal between left- and right-going half-scans, so separate 
templates are calculated for each case. The scan-synchronous 
contamination removed in this step is very small; Q and U 
maps typically change by 100^00 nK at the largest scales 
after its removal. 

The temperature T at each map pixel p is recovered from 
the filtered sum timestreams with 



(5) 



where w"*" is the weight assigned to each pair sum, and we 
have assumed that the effects of polarization leakage are neg- 
ligible. In other words, the temperature is obtained simply 
by binning the filtered detector timestreams into map pixels. 
Stokes Q and JJ are calculated from linear combinations of 
the difference timestreams using the matrix equation 




No.ofpairs ^ _ 

yep ^ 



No.ofpairs / ? n 



e(P) 
f/(p) 



(6) 



Here, w is the weight assigned to each pair difference, and 
a,/3 are PSB pair orientation angle factors defined as 

OLij = Hi cos 2ip2i,j - 721+1 COS 2ip2i+\j (7) 
Av = 72/sin2-02/j--72;+i sin2V'2/+ij-- (8) 

The 2x2 matrix in Equation[6]is singular for a single pair of 
V'2/,> and '02i+i,;, and the equation can be solved only by ac- 
cumulating more than one timestream sample in a given map 
pixel p. As p is observed with many detector polarization 
angles V', the off-diagonal a,;/3,7 terms average to zero, and 
the matrix becomes invertible. Although only two different 
polarization angles are required to invert the matrix, some in- 
strumental systematics average down as the number of obser- 
vation angles increases. By examining the determinant of the 
matrix, pixels (typically at the edge of the map) with insuffi- 
cient polarization angle coverage or low integration time are 
identified and masked from analysis. 

We choose the pair sum and difference weights 
to be proportional to the inverse variance of the filtered 
timestreams. The weights are evaluated from power spec- 
tral densities averaged over each set of azimuth scans (ev- 
ery 50 minutes), a period during which the noise proper- 
ties are approximately stationary. For each channel pair, the 
sum/difference weight for a scan set is calculated from the in- 
verse of the average value of the autocorrelation between 0.5 
and 1 Hz. 

5. MAP RESULTS 

Figure [2] shows the BiCEP maps of CMB temperature and 
Stokes Q and U parameters. The pixelization is 0.23°, cor- 
responding to Healpix nside=256, and the total observa- 
tion area covers ^ 2% of the sky. In the central part of the 



maps, the integration time exceeds 6x10^ detector-seconds 
per square degree, and the scan strategy smoothly apodizes 
the outer edges of the observed region. The temperature 
anisotropics are measured with high S/N and agree well be- 
tween the two frequencies. The rms noise per square de- 
gree, obtained from simulations, is 0.81 /iK and 0.64 /iK for 
100 GHz and 150 GHz, respectively, in the central region of 
the Q and U maps. 

Data from both frequency bands are combined to form the 
temperature, E, and B maps shown in Figure |3] The maps 
from each frequency band are smoothed to 1° resolution, and 
the left column shows the sum. Frequency jackknife maps 
that are formed by differencing the 100 and 150 GHz data are 
shown in the right column. The high S/N of the temperature 
measurement is demonstrated by the lack of structure in the 
frequency jackknife map. The faint striping, which is caused 
by residual unpolarized atmospheric noise, is successfully re- 
moved by PSB differencing. To form the E and B maps, the 
af,jj and a*,, coefficients are computed from apodized Q and 
U maps using the anafast utility in the Healpix code 
package. The coefficients are then boxcar filtered in a mul- 
tipole range of 70 < ^ < 280. Using the synfast utility, 
a^,jf are interpreted as Yi,„ coefficients from which the E and 
B maps are generated. The E frequency-sum map shows re- 
solved degree-scale structure of the expected amplitude, while 
the E and B jackknife maps, as well as the B signal map, are 
consistent with noise. These E and B maps are used only as a 
qualitative illustration of the S/N ratio of the data. Since the 
method of creating the maps has several shortcomings, e.g. 
there is some EIB mixing due to the finite survey area, we do 
not use the E and B maps in any quantitative analysis of the 
power spectra. 

6. POWER SPECTRUM ESTIMATION 

BiCEP employs two different analysis pipelines for power 
spectrum estimation: the first performs a curved-sky analysis 



using the publicly available Spice ( Chon et al.|2004] l pack- 
age to estimate the power spectra, and the second is a flat-sky 
pseudo-Q esti mator adapted from the QUaD pipeline that 
was used in the |Pryke et al.| (2009 1 main analysis. In the fol- 
lowing, we describe in detail the Spice pipeline. 

The power spectra are derived from inverse variance 
weighted T, Q, and U maps. Temperature and polarization 
correlation functions are computed from fast spherical har- 
monic transforms of the heuristically weighted maps. The 
polarization correlation functions are then decoupled so that 
EIB mixing that is caused by the finite survey area is removed 
in the mean. The correlation functions are apodized in or- 
der to reduce correlations in multipole space that result from 
incomplete sky coverage. Finally, estimates of the full sky 
power spectra are computed from integral transforms of the 
apodized correlation functions. The observed power spectra, 
as computed by Spice, are approximated as 



(9) 



where the superscripts, X, correspond to the six temperature- 
polarization combinations {TT, EE, BB, TE, TB, EB}. The 
convolution kernel, kJ^,, is the normalized Spice window 
function and depends on the apodization function applied to 
the correlation functions. The effect of timestream filtering 
is described by Ff, the ^-space transfer function; Bj = B^Hj 
is the product of the beam and pixel window functions, and 
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Fig. 3. — Data from BiCEP's 100-GHz and 150-GHz channels are combined to form temperature, E, and B signal and jackknife maps. The E and B maps 
are apodized to downweight noise-dominated edge pixels. The temperature anisotropics are measured with high S/N, and the E signal map shows resolved 
degree-scale structure. The B signal map and the E and B jackknife maps are consistent with noise. 



Nf is the power spectrum of the noise convolved with the 
Spice kernel. In the following subsections, we describe the 
steps in solving Equation [9] and recovering an estimate of the 
underlying power spectrum, Cf . 

6. 1 . Noise Subtraction 

The first step in power spectrum estimation is calculating 
and subtracting the noise bias N^. We estimate with Monte 
Carlo simulations of instrument noise: starting from a noise 
model, we generate simulated timestreams that are filtered, 
co-added into maps, and processed with Spice. The result- 
ing simulated noise spectra are averaged over many realiza- 
tions to form A'^. 

The BiCEP noise model is derived from gain-adjusted PSB 
pair sum and difference timestreams (Equation Hh under the 
assumption that the timestream S/N is negligible, i.e., the 
signal is the noise. The S/N is < 10% and < 0.1% for the 
sum and difference timestreams, respectively. The sum and 
difference timestreams are polynomial-filtered and Fourier 
transformed over each half-scan, and all auto-correlations 
and cross-correlations between channel pairs are computed to 
form the complex frequency-domain noise covariance matrix 



N(/). To form the noise model, we average N(/) over the 
100 half-scans within each scan set and then average into 12 
logarithmically spaced bins spanning 0.05-5 Hz. 

To construct simulated correlated noise timestreams, we 
take the Cholesky decomposition N(/) = L(/)L^^(/) of the 
noise covariance matrix and multiply a vector of pseudo- 
random complex numbers p{f) by L(/): 



v(/) = L(/)p(/). 



(10) 



The complex numbers in p{f) have Gaussian-distributed real 
and imaginary components and are normalized such that the 
magnitude has a standard deviation of one. The resulting 
product v(/) has the same statistical properties as the data 
and is inverse Fourier transformed to obtain a set of simu- 
lated noise timestreams. Scan-synchronous templates are cal- 
culated and subtracted from each set of azimuth scans, and the 
filtered noise timestreams are then co-added into maps. 

The noise bias, A^^, is estimated by averaging the power 
spectra from an ensemble of noise-only maps. Figure|4]shows 
(red dashed lines), calculated from the average over 500 re- 
alizations, in comparison to raw power spectra. BiCEP mea- 
sures TT and TE with high S/N, and the noise contribution 
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is negligible up to ^ 330. The noise is mostly uncorrelated 
between temperature and polarization, so for TB is dis- 
tributed around zero, and the same is true for EB. In con- 
trast, noise comprises the bulk of the EE and all of the BB 
amplitude at ^ > 150, which illustrates the need for careful 
noise modeling and subtraction. We have studied the impact 
of n oise misestimation by intentionally scaling the noise spec- 
tra ( Takahashi et al.|2010 1, and we find that the largest possi- 
ble error in A^^ is ±3%. 

6.2. Beam and Pixelization Corrections 

The finite resolution of the telescope and pixelization of 
the maps result in suppression of the observed power spec- 
tra at small angular scales. This suppression of power is de- 
scribed by B\ = B\H^, the product of the beam and pixel win- 
dow functions, and is illustrated in FigureJS] In this analysis, 
we approximate Bg as the Legendre transform of the average 
beam within a given frequency band, which is assumed to be 
a circular Gaussian. The full widths at 100 and 150 GHz are 
0.93° and 0.60°, respectively. The pixel window function 
is supplied by Healpix and corresponds to nside=256. In 
the case of BiCEP, Bi varies slowly with respect to the Spice 
kernel, kJ^, , and can therefore be pulled out of the convolu- 
tion. The observed power spectra, after noise subtraction, are 
divided by B\ to correct for the effects of beam suppression 
and pixelization: 



{Cj-Nf)/B\c^Y.^, 



X T7X /^X 



(11) 



and denote the first and second derivatives in spherical 
coordinates. Assuming Gaussian beams for each frequency 
band. Equation 12 reduces to 



c/(p) = m(p) + Vm(p)(ri - p) + 
1 



rTr 



(13) 



where Ai/i and A6' are the longitude and latitude offsets be- 
tween the pointing vector and the pixel center p. We apply 
Equation [13] to simulate signal-only detector timestreams ac- 
cording to BiCEP's scan strategy. Measured PSB pair cen- 
troids, detector orientation angles, and cross-polar leakage 
values are included in the simulations. For the purpose of 
characterizing the effects of timestream filtering, all differen- 
tial beam systematic effects are turned off so that there is no 
mixing between temperature and polarization. The simulated 
timestreams are filtered and weighted in exactly the same way 
as the real data and then co-added into maps. 

Once we have "BiCEP-observed" signal-only maps in hand, 
we compute the power spectra with Spice and average the 
results over many realizations. Because the input CAMB spec- 
tra have r = 0, any non-zero BB power in the simulation out- 
puts is interpreted as contamination from E-vaodt polarization 
that is induced by timestream filtering. We therefore apply a 
correction. 



(c; 



BB \tBB /^j 

lyf — u, 



BB\ 



: only)/ BI 



EBB pBBf^BB 



(14) 



6.3. Filter Corrections 

The bolometer timestreams are cleaned by subtracting a 
third-order polynomial and scan-synchronous template; this 
cleaning procedure also has the effect of removing large-scale 
CMB signal. The amount of signal loss is described by Fg, 
the £-space transfer function imposed by timestream filtering. 
In addition, although the Spice estimator is unbiased in the 
mean, the timestream processing removes spatial modes from 
the observed T, Q, and U maps, generically introducing cou- 
plings between the observed E and B spectra. These couplings 
are small, but are of interest for the leakage of the relatively 
large /s-mode signal into B-mode polarization. We character- 
ize both F( and E-to-B leakage with signal-only simulations. 

The signal simulation procedure begins by generating 



model power spectra using the CAMB (Lewis et al. 
software package. We use ACDM pararneters derrvei 



from 



Wmap five-year data ( [Hinshaw et al.|[2009| l and a tensor- 
to-scalar ratio of zero. From the model spectra, we use the 
synf ast utility in Healpix to create an ensemble of sim- 
ulated CMB skies pixelized at 0.11° resolution (nside=512) 
and convolved with the BiCEP beams. Actual pointing data 
are used to generate smoothly interpolated PSB timestreams 
from the simulated T, Q, U maps and their derivatives. A PSB 
timestream sample ii(p) that falls into a pixel centered at p 
is expressed as a convolution of the beam 7'(r-r/,), which is 
centered at r/,, with a second-order Taylor exp ansion of the 
sky signal m around p ( Hivon & Ponthieu|2010[ l: 



d{p)= / drV(r-ri,) m(p) + Vm(p)(r-p) 



1 



+-(r-p)'DMp)(r-p) 



(12) 



to the BB power spectrum of the data, where C^^Iee only is 
the ensemble average BB spectrum from the r = signal sim- 
ulation outputs. The amplitude of this correction is roughly 
£{£+ l)Cf\EE oniy/(27r) ~ 3 X 10-3 at £ - 100, compa- 
rable to inflationary BB power for r = 0.05, and the sample 
variance is about 3 x lO""^ fiK^. The correction factors for the 
other spectra are negligibly small, so only the BB spectrum is 
adjusted with this procedure. 

To quantify the suppression of large-scale power from 
timestream filtering, we begin by assuming that Ff, like Bj, 
varies slowly in comparison to the Spice kernel and can be 
pulled out of the convolution. We choose to calculate the 
binned filter suppression factor, F^, defined as 



'-F,^^Pl 



2tt ■ 



^KU'Cf,, 

(15) 

where the binning operator is a top hat. (In the case of the 
BB spectrum, the left hand side of the equation also contains 
C^^Iee only, as in Equation 14 ) We obtain the filter suppres- 



Here, m(p) = g[r(p) + 7(2(p)cos2-i/' + t/(p)sin2V')], and V 



sion factors by comparing me power spectra of the BlCEP- 
observed signal-only maps to those of the input synfast 
maps; F^ is the ratio of the spectra, after multiplying by 
£(£+ l)/(27r) and binning. 

FigureBlshows F^ averaged over 500 signal simulations us- 
ing BiCEP s timestream filtering. In each simulation, signal- 
only timestreams are generated from the full two years of 
pointing data. At ^ ^ 100, the value of F^ is about 0.7 for 
all spectra and rises slowly as £ increases. Identical filter 
functions are used for the EE and BB spectra, and the filter 
functions of the cross-spectra are calculated as the geometric 
mean of those determined from the auto-spectra. The validity 
of this approach for obtaining the BB and cross-spectra filter 
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Fig. 4. — Raw power spectra, uncorrected for noise and filter bias, are shown liere for 500 signal-plus-noise simulations (gray lines) and BiCEP 150 GHz data 
(black points). The dashed red lines indicate the level of noise bias, £(£+ \ )Ni/(2tt), and the black lines show the ACDM spectra used as inputs for the signal 
simulations. Noise dominates much of the EE spectrum and all of the BB spectrum. The effect of filter bias is visible in the low-£ portion of the TT spectrum, 
where the simulated spectra fall below the model curve. 
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Fig. 5. — Signal-only simulations are used to evaluate the filter function 
, shown here for 150 GHz {Fj^ for other frequency combinations look 
similar). The beam-pixel window functions S| = B^^H^ are shown for 100 
and 150 GHz. 



The uncertainties in the power spectra consist of two com- 
ponents, one that is proportional to the signal itself (sample 
variance), and another that depends on the instrumental noise. 
We estimate the errors by examining the variance of power 
spectra from simulated signal-plus-noise maps, which exactly 
encode time-dependent correlated noise, scan strategy, and 
sky coverage. We add t he si mula ted noise-only and signal- 
only maps, described in ^6.1| and ^6.3[ and power spectra are 
calculated for each realization using the same A^^, B\, and 
as applied to the real data. If the simulations include a reason- 
able model of the signal and faithfully reproduce all the noise 
properties of the experiment, then the data and simulations 
should be indistinguishable. Figure |4] shows raw power spec- 
tra, uncorrected for noise and filter bias, from 500 signal-plus- 
noise simulations at 150 GHz. The raw spectra of the actual 
data, shown by the black points, lie within the scatter of the 
simulations. We calculate the band power covariance matrix 
from the ensemble of simulations, after applying noise, filter, 
and beam corrections. The band power errors are obtained 
from square root of the diagonal terms of the matrix. 



functions has been confirmed with simulations over a limited 
multipole range. Dividing by is the final step in obtaining 
the BiCEP band power estimates: 



J{l+\){Cj-Nf) 



In 



c^BB _ \ ' pb 



,Xj^BB (16) 
"""■^ (17) 



ea+i) (cf-jvf-crUfi only) 

6.4. Error Bars 



6.5. Band Power Window Functions 

In order to compare our band power estimates to a theo- 
retical model, we need a method to calculate expected band 
power values from theoretical power spectra. The relationship 
between the model and the expected band powers is described 
by band power window functions, ujf , defined as 



(18) 



10 



CHIANG ET AL. 




300 350 
Multipole / 



F ig. 6 . — BiCEP's band power window functions are defined by Equa- 
tion 1 8 and this figure shows u)'l/l. 



where "^i = £(£+ l)Q/(27r). The window functions are given 
by 



2£ 



and are normalized such that 



E 



(19) 



(20) 



The window functions depend on the Spice kernel, which 
depends on the apodization function applied to the correla- 
tion functions. For this analysis, we apodize the correlation 
functions with a cosine window that spans 50°. We choose 
to use a uniform binning of width = 35, spanning a multi- 
pole range of 21 < £ < 335. This choice of bin width provides 
minimally correlated band powers while preserving the spec- 
tral resolution determined by the width of the Spice kernel. 
Figure|6]illustrates BiCEP's band power window functions for 
the nine I bins. 

Because the BB power spectrum is debiased with the proce- 
dure described in Equation [14] we set the EE-to-BB window 
functions to zero. This is a valid approximation as long as the 
EE spectrum of the signal under consideration is statistically 
consistent with the measured BiCEP band powers (which, as 
shown in 5 12 are well described by a concordance ACDM 
cosmology )7'^ 

7. POWER SPECTRUM RESULTS 

Figure |7] shows the full set of BiCEP spectra plotted with a 
ACDM model derived from Wmap five-yeai" data. The 100 
and 150-GHz auto-spectra are shown, as well as the lOOx 150 
cross-spectra. In the case of the TE, TB, and EB spectra, we 
also show the 150x 100 cross-correlation. For each spectrum, 
we present nine band powers with a uniform bin width of A£ = 
35, spanning 21 < i < 335. The TT, TE, and EE spectra 
are detected with high significance and are already sample- 
variance limited, and there is no detection of signal in BB, 
TB, and EB. The results from BiCEP's two analysis pipelines 



agree well with each other, and Figure[TT|shows a comparison 
of the frequency-combined power spectra (^TT). 

As a cross-check, we have also derived TT, TE, and TB 
spectra using Wmap five-year temperature data in BiCEP's 
CMB field (open circles in Figure [7]). The Wmap temper- 
ature maps are smoothed and filtered identically to BiCEP, 
as described in 5 3.2 The TT points are calculated from the 
cross-correlation of the WMAP Q and V-band maps, and the 
TE and TB spectra are calculated using Wmap V-band tem- 
perature data and BiCEP polarization data. For this compar- 
ison, we do not subtract noise bias and instead rely on the 
fact that the pairs of maps have uncorrelated noise. We also 
do not attempt to assign error bars. Qualitatively, the spec- 
tra formed using Wmap temperature data agree well with the 
spectra from BiCEP temperature and polarization data. Both 
the BiCEP and Wmap temperature maps are strongly signal- 
dominated; apparently the differences between them, includ- 
ing the noise as well as potential systematics, are at a level that 
has little impact on the TT, TE, or TB power spectra results. 

8. CONSISTENCY TESTS 

8.1. Jackknife Descriptions 

We check the self-consistency of the power spectra by per- 
forming jackknives, statistical tests in which the data are split 
in two halves and differenced. The split is performed at the 
mapmaking stage, and the resulting differenced map should 
have power spectra that are consistent with the expected resid- 
ual signal level (nearly zero) after subtracting noise bias. The 
interaction of timestream filtering with the details of the split 
causes imperfect signal cancellation when forming jackknife 
maps, but in practice, this residual signal is small. 

The data are tested with six jackknives that are sensitive to 
different aspects of the instrument's performance. In the scan 
direction jackknife, the data are split into left- and right-going 
azimuth half-scans. Failures generally point to a problem in 
the detector transfer function deconvolution, or thermal insta- 
bilities created at the scan endpoints. The elevation coverage 
jackknife is formed from the two CMB observations in each 
48-hr cycle; each observation covers the same azimuth range 
but starts from a different elevation. This jackknife is sensi- 
tive to ground-fixed or scan-synchronous contamination. Bi- 
CEP observes at four fixed boresight orientation angles that 
can be split into two pairs, {-45°, 0°} and {135°, 180°}, to 
form a boresight angle pair jackknife. This test is perhaps the 
most powerful of the jackknives performed and is sensitive to 
many factors, including thermal stability, atmospheric opac- 
ity, relative gain mismatches, differential beam pointing, and 
ground pickup. In the temporal jackknife, the 8-day observing 
cycles — 48 hr at each of the four boresight angles — are inter- 
leaved to form the two halves, and this jackknife tests sensi- 
tivity to weather changes. The season split jackknife simply 
divides the data into the two observing seasons, and failures 
reflect any changes made to the instrument between the two 
years. In particular, the focal plane thermal architecture was 
improved for the 2007 season, and the temperature control 
scheme was changed. The focal plane QU jackknife splits 
the detectors into two groups according to their polarization 
orientation within the focal plane (approximately alternating 
hextants) and is a method of probing instrumental polarization 
effects. 

Power spectra of jackknife maps are computed with the 
method described in ^ using simulated jackknife noise and 
signal-plus-noise maps to subtract noise bias and assign er- 
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Fig. 7. — Black points show the full set of BiCEP's power spectra. The horizontal axis is multipole moment I, and the vertical axis is l)Q/(27r) in units 
of /iK^. The spectra agree well with a ACDM model (black lines) derived from Wmap five-year data and r = 0. The gray points coirespond to the boresight angle 
pair jackknife; note that although the TT , TE, and TB jackknife failures are statistically significant, the amplitudes are small compared to the signal. The open 
circles show spectra calculated using WMAP five-year temperature data, smoothed and filtered identically to BiCEP. We have cross-correlated the WMAP Q and 
V-band data in BiCEP's CMB field to obtain the TT points, and the TE and TB spectra are calculated from the cross-con'elation of WMAP V-band temperature 
data with BiCEP polarization data. 



12 



CHIANG ET AL. 



ror bars. The filter function is applied so that the mag- 
nitude of any non-zero jackknife band powers can be com- 
pared to the amplitude of the non-jackknife spectra. For each 
jackknife spectrum, we calculate over nine bins spanning 
21 < £ < 335. To account for the expected level of residual 
signal, the values are evaluated with respect to the aver- 
age jackknife spectra from an ensemble of signal-only simu- 
lations. The criteria for jackknife success or failure are based 
on the probability to exceed (PTE) the value, which is cal- 
culated from the distribution of in 500 signal-plus-noise 
simulations. Jackknife victory is declared when (1) none of 
the PTEs are abnormally high or low, given the number of x^ 
tests performed, and (2) the PTEs are consistent with a uni- 
form distribution between zero and one. 

8.2. Jackknife Results 

The jackknife spectra for the six different data splits appear 
similar in that the band powers are all distributed around zero. 
We therefore show only the spectra from the boresight pair 
jackknife (Figure [7] gray points), which is arguably the most 
stringent of the six tests. Table[T]lists the PTE values from all 
of the jackknife x^ tests for BiCEP's polarization-only spectra 
(EE.BB.EB). Out of all the PTEs, the lowest value is 0.014, 
in the 150 GHz focal plane QU jackknife. This low value 
is, however, consistent with expectations from uniformly dis- 
tributed PTEs over the 60 polarization-only jackknives. Fig- 
ure |8] shows that the PTEs are consistent with a uniform dis- 
tribution between zero and one. 

The polarization jackknife tests are the most powerful probe 
of the accuracy of the noise model. In addition to the x^ tests, 
we have also expressed the jackknife spectra as band power 
deviations with respect to the mean of signal-only simula- 
tions. The sum of the band power deviations in each jackknife 
spectrum provides an additional and more precise gauge of 
the correct estimation of instrument noise. We have verified 
that the band power deviation sums in the data are consistent 
with signal-plus-noise simulations and are not systematically 
biased high or low, thus confirming that the noise levels are 
correctly estimated. Furthermore, we have probed the sensi- 
tivity to systematic variations in the modeled noise amplitude 
over a range that is comparable to the S/N of the sum and 
difference timestreams, and find no significant changes in the 
jackknife spectra. The relative insensitivity of the jackknifes 
to the amplitude of the noise model validates the procedure 
described in §6.1| 

In contrast to the polarization data, the temperature data dis- 
play significant jackknife failures. There is an excess of small 
PTE values in the TE and jackknives, and most of the TT 
jackknife PTEs are smaller than 0.002, which is the resolu- 
tion from 500 simulations. The TT, TE, and TB jackknife 
PTEs are therefore not listed in Table [T] We attribute these 
failures to the fact that BiCEP's temperature maps have high 
S/N (see TT plot in FigureHll, and the jackknives are therefore 
extremely sensitive to small gain calibration errors or imper- 
fections in modeling and subtracting unpolarized atmospheric 
emission. (As described in ^6.1| the < 10% S/N in the pair- 
sum timestreams is a known imperfection of the BiCEP noise 
model for temperature data.) Although the TT , TE, and TB 
jackknife failures are statistically significant. Figure [7] illus- 
trates that the amplitudes of the jackknife spectra are small 
compared to both the amplitude and errors of the signal spec- 
tra. The magnitudes of the TT and TEITB jackknife band 
powers are typically 1-10 jiK? and 0.1-1 jjX?, respectively. 



TABLE 1 

Jackknife PTE Values from Tests 



Jackknife 


100 GHz 


150 GHz 


100x150 


150x100 


Scan direction 
EE 
BB 
EB 


0.532 
0.640 
0.816 


0.588 
0.568 
0.962 


0.740 
0.212 
0.924 


0.358 


Elevation coverage 
EE 
BB 
EB 


0.576 
0.584 
0.872 


0.546 
0.288 
0.728 


0.924 
0.618 
0.892 


0.892 


Boresight angle 
EE 
BB 
EB 


0.916 
0.242 
0.912 


0.448 
0.548 
0.100 


0.320 
0.592 
0.392 


0.944 


Temporal split 
EE 
BB 
EB 


0.378 
0.788 
0.370 


0.208 
0.020 
0.580 


0.796 
0.852 
0.476 


0.232 


Season split 
EE 
BB 
EB 


0.564 
0.790 
0.806 


0.716 
0.992 
0.514 


0.216 
0.056 
0.456 


0.986 


Focal plane QU 
EE 
BB 
EB 


0.670 
0.896 
0.236 


0.014 
0.804 
0.806 


0.994 
0.576 
0.234 


0.560 




0.6 0.7 0.8 0.9 1 
EE, BB, EB jackknife PTE 

Fig. 8. — Probabilities to exceed the values from EE, BB, and EB 
jackknife tests are consistent with a uniform distribution between zero and 
one. 

In all cases, the error bars of the non-jackknife spectra are at 
least a factor of 10 larger, except in the highest t bin. 

We have performed the same jackknives with both analysis 
pipelines, and the results are in excellent agreement. The al- 
ternate pipeline confirms that the polarization data pass the x^ 
tests, and although the TT, TE, and TB spectra do not pass, 
the amplitudes of the jackknife spectra are small. 

9. SYSTEMATIC UNCERTAINTIES 

Systematic errors that arise from uncertainties in instrument 
characterization can be separated into two categories: (1) er- 
rors that mix temperature, fi-mode, and Z?-mode polarization; 
and (2) calibration uncertainties that affect only the scaling or 
amplitude of the spectra. In this section, we describe the dom- 
inant sources of systematic error in BiCEP and the expected 
impact on the power spectra; these systematics are summa- 
rized in Figure |9] A complete description of all potential Bl- 
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Fig. 9. — Dominant systematic uncertainties in BiCEP's band powers are 
compared to tlie total statistical eiTor, which comes from sample variance 
and noise. The horizontal axis is multipole moment £, and the vertical axis 
is £(i + 1 )Q /(2-k) uncertainty in units of fiK^. The contributions from the 
signal-noise cross terms to the total statistical error in TE are not shown. 
We have assumed that gain mismatch and differential pointing systematics 
leak temperature to EE and BB in equal amounts. 



CEP systematics and the methodology for propagating the er- 
rors to the power spectra is gi ven in the accompany ing instru- 
ment characterization paper ( [Takahashi et al.|2010| l. 

9. 1. Temperature and Polarization Mixing 

We are mostly concerned with systematic errors that mix 
the bright temperature signal into polarization and thus in- 
duce a false B-mode signal. We define a benchmark for each 
systematic error such that the false B-mode amplitude is no 
greater than the peak of the inflationary BB spectrum. For 
BiCEP's target of r = 0.1, this requirement corresponds to 



£(£+ l)Cf/(2n) = 7 X lO'^ /iK^ at . 



100. 



The primary source of potential temperature leakage into 
polarization is differential gains within PSB pairs. Gain mis- 
match effects can also arise from other systematics, such 
as errors in the bolometer transfer functions, which act 
as frequency-dependent gains. We have characterized the 
common-mode rejection of PSBs by examining the cross- 
correlation of pair-sum and pair-difference maps. There is 
no statistically significant evidence for gain mismatch in the 
data, and we set an upper limit of A{gi /g2)/{gi/g2) < 1-1% 
on differential gains. 

A second source of potential temperature and polarization 
mixing is beam mismatch within PSB pairs. We describe 
beam mismatch with three quantities: differential beam size, 
pointing, and elUpticity. Of these three, the dominant ef- 



fect in BiCEP is differential pointing, which is stable over 
time and has been measured with a median amplitude of 
(ri -r2)/a = 1.3 ±0.4%, where a is the average Gaussian 
beam size within a pair of PSBs. Measured upper limits on 
differential size and ellipticity are negligible. 

Most systematic errors interact with the scan strategy in 
complex ways, and the exact effects on the power spectra can 
be computed only thro ugh signal sim ulations. We fo llow the 
formalism p resented in Hivon & Ponthieu (20l0| and Shimon 
|et al.| (|2008) to calculate the expected level of lalse BB in Bl- 
CEP. Starting from synf ast maps with r = 0, Equation 13 
is used to generate simulated signal timestreams that inclime 
the effects of gain mismatch and differential pointing. The 
timestreams are filtered and co-added into maps, and the am- 
plitude of the BB power spectrum at ^ ~' 100 is compared with 
the 7 X lO"-' fiK^ benchmark. 

In the differential gain simulations, we randomly assign 
1.1% rms gain mismatch to PSB pairs, and we find that the 
expected false BB amplitude is 7.4 x lO"-' /zK^ and 2.9 x 
10-2 fiK^ at £ - 100 for 100 and 150 GHz, respectively. Al- 
though this amplitude exceeds the r = 0. 1 benchmark. Fig- 
ure l9]illustrates that it is small compared to the statistical error 
in tms analysis of the two-year data. In a future analysis of 
the entire BiCEP data set, we anticipate placing tighter con- 
straints on PSB gain mismatch as the noise levels integrate 
down. With further work, we are confident that uncertainty in 
gain mismatch can be substantially reduced from the current 
1.1%. 

Differential pointing has been precisely characterized for 
each of BiCEP's PSB pairs, so we run signal simulations using 
the measured centroid offset vectors, rather than randomized 
distributions. The false BB from differential pointing has an 
e - 100 ampUtude of 2.7 x 10"^ ^iK^ and 4.2 x 10"^ fiK^ at 
100 and 150 GHz, respectively. These amplitudes are slightly 
smaller than the r = 0.1 benchmark and are well below the 
noise level of the initial two-year data set. In a future analy- 
sis, it may be possible to use the measured centroid offsets to 
correct for systematic effects. 

We emphasize that, in addition to the differential gain 
and pointing discussed here, most uncertainties in instrument 
characterization create false positive BB signal. The fact that 
BiCEP's BB spectra are consistent with zero and pass jack- 
knives demonstrates that we have achieved sufficient control 
over systematic errors in this analysis. Furthermore, until a 
positive B-mode detection is made, the presence of systematic 
effects that produce spurious polarization could only make the 
reported BB upper limits higher (more conservative) than they 
would be otherwise. 

9.2. Absolute Gain and Beam Uncertainty 

The scaling of the power spectra is determined by the ab- 
solute gain factors that convert detector uni ts to tempera- 
ture, and the 2% uncertainty in this gain ( ^3.2[ ) translates 
into a 4% uncertainty in the power spectrum amplitude. The 
polarized spectra have additional amplitude uncertainty that 
arises from errors in the cross-polar leakage. A systematic 
error of Ae = 0.01 corresponds to 3.9% amplitude uncer- 
tainty for polarization-only spectra {EE, BB, EB) and 2.0% 
for temperature-polarization cross-spectra (TE, TB). There- 
fore, the combined amplitude uncertainty = iSC^jCf 
from gain and e errors is 4% for TT , 4.5% for TEITB, and 
5.6% for EEIBBIEB. These uncertainties are shown in Fig- 
ure |9] assuming ACDM spectra with r = 0. 
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For Gaussian beams, measurement errors in the beam 
widths introduce a fractional uncertainty, 



Si EE AQ /Cf = e 



■X if^X _ cr^£(^+l)(5-+2(5)_ 



1-5, 



(21) 



in the power spectrum amplitude as a function of I. Here, 
a = FWHM/ v^8Tn(2), and S is the fractional beam width er- 
ror Aa/a. Because this band power uncertainty is degenerate 
with absolute gain error, we subtract the mean, S, calculated 
over 56 < ^ < 265, the angular scales over which we perform 
absolute calibration. For this analysis, we use average beam 
widths of 0.93° and 0.60° at 100 and 150 GHz, respectively. 
Although the distribution of beam widths in the focal plane 
varies by ±3%, the measurement precision is ±0.5%; we 
therefore expect the effective S to lie somewhere in between. 
We calculate the maximum beam width error allowed by our 
calibration cross-spectra (Equation |2]i, which are very close 
to flat, and we consti-ain S < 1.2% and 6 < 2.8% at 100 and 
150 GHz, respectively. Figure [9] shows the expected power 
spectrum errors from these maximum allowed beam uncer- 
tainties, which are most likely conservative. The systematic 
error is smaller than the statistical error in all cases except the 
TT spectra at high £, where the levels of statistical and beam 
systematic uncertainty are comparable. 

10. FOREGROUNDS 

The BiCEP CMB region was chosen to have the lowest fore- 
ground dust emission for a field of that size, and we do not 
expect foreground contamination from dust or other sources 
to be significant at the current depth in the maps. To verify 
this, we estimate the levels of contamination in BiCEP CMB 
data from three potential foreground sources: thermal dust 
emission, synchrotron radiation, and extragalactic radio point 
sources. 

10.1. Thermal Dust 

Polarized dust emission in the BiCEP CMB field is esti- 
mated from "EDS Model 8" (Finkbeiner et al. 1999). We as- 
sume a 5% polarized fraction, guided by a study of Wmap 
data that shows that hi gh-latitude dust has a mean fractional 
polarization of 3.6% (Kogut et al. 2007 1. The dust tem- 
perature and polarization model is extrapolated to 100 and 
150 GHz, and filtered according to BiCEP's scan strategy. 
The resulting polarized dust emission is t{i + \)C ^ / {Itt) = 
9.6 X 10-5 /iK^ and 6.1 x lO"'* at £ = 100 for 100 and 
150 GHz, respectively. These values are 2 orders of magni- 
tude below the 95% confidence level for upper limits on the 
BB amplitude from BiCEP (discussed in ^T3] l, as shown in 
Figure [To] 

We nave also tested for dust contamination in the BiCEP 
maps by studying the cross power spectrum between EDS 
Model 8 maps and BiCEP temperature and polarization maps. 
The cross power spectrum of the real data is consistent with 
the distribution from signal-plus-noise simulations, providing 
additional evidence that the BiCEP maps contain no spatial 
correlation with the EDS dust maps. 

10.2. Synchrotron 

To estimate the polarized synchrotron emission in the Bi- 
CEP field of view, we have used the Wmap Ma rkov Chain 
Monte Carlo (MCMC) synchrotron maps ( Gold et al.|2009| l, 
extrapolated to 100 and 150 GHz using the mean spectral iii- 
dex of the BiCEP field of view, calculated from the spectral 
index maps provided in the same analysis. As for the EDS 
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Fig. 10. — Expected levels of polarized dust and polarized synchrotron in 
the BiCEP CMB field, assuming FDS Model 8 with 5% polarization fraction 
and Wmap MCMC polarization maps extrapolated to 100 and 150 GHz. 
These estimat ed fo reground levels are much lower than the BB upper limits 
presented in ^13| Theoretical BB curves for r = 0.1, with and without 
gravitational lensmg, are also shown for comparison. 



dust maps, we filter the extrapolated 100 and 150 GHz syn- 
chrotron maps with the BiCEP scan strategy, and find the es- 
timated level of polarized synchrotron emission at £ ^ 100 
to be 3 X 10"^ AiK^ and 4 x 10"^ yuK^ for 100 and 150 GHz, 
respectively, both below the level of BiCEP sensitivity. Fur- 
thermore, the Wmap MCMC map has poor S/N in regions far 
from the Galactic plane, and we have found that within the 
BlCEP CMB field, the map is dominated by variance in the 
Monte Cai-lo fit (see Gold et al. ,200 9, for details). The esti- 
mated levels of polarized synchrotron emission in Figure [10] 
should therefore be viewed as conservative estimates or up- 
per limits. We have also derived synchrotron estimates from 
Wmap K-band polarization maps and the temperature maps 
of Haslam et al. ( 198 1 ), assuming a 30% polarization fraction; 
both methods give synchrotron estimates that are lower than 
tiiose from the Wmap MCMC maps. 

Similarly to our analysis for thermal dust, we have stud- 
ied the cross power spectrum of the synchrotron maps with 
the BiCEP maps and found there to be no significant spatial 
correlation of the BiCEP data with the synchrotron emission. 

10.3. Point Sources 

At degree-scale resolution, the BiCEP maps do not show 
any obvious point source detections, so we rely on a com- 
binati on of die 4.85-GHz Pai-kes-MlT-NRAO (PMN) sur- 
vey (Wright et al. 1994 1, the Wmap point source cata- 



log (Wright "^ al.||2b09[ ),"and the ACBAR catalog ( jReichardt) 
|et al.||2009) ) to search for point source contamination. We 
search for point source contamination by optimally filtering 
CMB and noise fluctuations out of the BlCEP temperature 
map and determine the significance of the resulting pixel val- 
ues by repeating the process with simulated maps of CMB 
with detector noise. Although the resulting maps have a few 
2cr detections at the suspected point source locations, there is 
no statistical evidence for point source contamination above 
the expected Gaussian distribution of noise. As a further test, 
we have simulated the effects of masking out the 27 ACBAR 
sources that lie within the BlCEP field and have found it has 
no significant impact on the power spectra. 
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10.4. Frequency Jackknife 

The CMB and foreground emission have different fre- 
quency dependence, so we can test for the presence of fore- 
ground contamination in BiCEP data by performing a fre- 
quency jackknife. We difference the 100 and 150-GHz maps, 
compute the power spectra, calculate and compare the re- 



sults to signal -plus-noise simulations, as described in 5 8.1 
The probabilities to exceed the values are {0.050, 0.152, 
0.732} for EE, BB, and EB, respectively. We find no evidence 
for foreground contamination in the frequency jackknives. 

11. COMBINED SPECTRA 

We combine the spectra from the different observing fre- 
quencies by taking a weighted average for each band power. 
To obtain the weights, we use signal-plus-noise simulations to 
calculate the covariance matrices from the various frequency 
combinations (100, 150, 100x150, and 150x100 GHz). 
There are three unique combinations for TT, EE, and BB, 
and four combinations for the other spectra. The weights are 
calculated from the row/column sums of the inverse of the 
covariance matrices. The error bars of the combined spec- 
tra are determined by applying the same combination weights 
to signal-plus-noise simulations. For fully noise-dominated 
spectra (such as BB), the error bars of the combined spectra 
improve by IQ%-4Q% compared to the errors from 150 GHz 
data alone. 



As suggested by Bond et al. ( 2000 1, we apply a transforma- 
tion 

Zh = \n(%+xi,) (22) 

to account for the fact that the probability of the true model 
value, given an observed band power, is offset-lognormally 
distributed. The offsets x/, describe the noise power spectra 
on the sky (i.e., corrected for filter and beam bias) and are 
calculated from simulations. We calculate Xh for the TT, EE, 
and BB spectra, but we assume Gaussian distributions for the 
TE, TB, and EB band powers since the values can be nega- 
tive. The BiCEP band powers, Xh offsets, covariance matri- 
ces, and band power window functions are available online at 
http://bicep.caltech.edu 

Figure [TT] shows a comparison of the frequency-combined 
spectra \vltn a ACDM model derived from Wmap five-year 
data. The power spectrum results are confirmed by the alter- 
nate analysis pipeline (open circles. Figure [TT). BiCEP con- 
tributes the first high S/N polarization measurements around 
i ~ 100, as illustrated by Figure 12 which shows the EE peak 
at £ ~ 140 in greater detail; the BB spectrum is overplotted for 
comparison. Figure 13 shows BiCEP's TE and EE spectra, as 
well as the 95% confidence upper limits on BB, in addition to 
other recent CMB polarization data. To obtain the BB upper 
limits, we apply offset-lognormal transformations to the band 
powers and integrate the positive portion of the band power 
probability distributions up to the 95% point. BiCEP mea- 
sures EE in a multipole window that complements existing 
data from other experiments, and all nine band powers have 
> 2(7 significance. The constraints on BB are the most pow- 
erful to date. 

12. CONSISTENCY WITH ACDM 

The power spectra of the CMB are well described by a 
ACDM model, which, in its simplest form, has six parame- 
ters that have been constrained by numerous experiments. We 
check the consistency of the BiCEP band powers with this 
model by performing a test. We start by using CAMB to 



calculate theoretical power spectra, using ACDM parameters 
derived from Wmap five-year data (and r = 0), and we then 
compute expected band power value s, "j^ f , using the band 
power window functions described in ^6.5| Absolute gain and 
beam systematic errors (G^ and S/,, as described in 5 9.2 1 are 



included by adding their contributions to the band power co- 



variance matrix, M^^: 



M 



: Mf„ + (G" )^"^/ + SAK %. 



(23) 



The Sh factors are formed from linear combinations of the 
four frequencies (100 GHz auto, 150 GHz auto, 100x150, 
150x 100), using the weights described in ^TT] Because M^^ 
is obtained from a limited number of simulations, the far off- 
diagonal terms are dominated by noise; we therefore use only 
the main and first two off-diagonal terms of M^^ in this calcu- 
lation. (We have tested that results are essentially unchanged 
including one, two, or all off-diagonal terms.) For each power 
spectrum, the observed and theoretical band powers are com- 
pai^ed by evaluating 



X 



(24) 



over the nine bins that span 21 < ^ < 335. In the case of the 
TT , EE, and BB spectra, offset-lognormal transformations. 



Zf = lnC^/+4), 
Z^ = ln(<+4), 



(25) 
(26) 
(27) 



are applied to the data, expected band powers, and inverse 
covariance matrix, and is calculated using the transformed 
quantities. 

We perform the same calculations for a set of 500 signal- 
plus-noise simulations, and the simulated x^ distributions are 
used to determine the probabilities to exceed the x^ values 
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of the data. The x^ and PTE values are listed in Figure 
which shows a comparison of our data with the ACDM model 
The BiCEP data are consistent with ACDM, and this result is 
confirmed by the alternate analysis pipeline. 

13. CONSTRAINT ON TENSOR-TO-SCALAR RATIO FROM BB 

BiCEP was designed with the goal of measuring the BB 
spectrum at degree angular scales in order to constrain the 
tensor-to-scalar ratio r. We define r = A^(fc()) / A^(fc()), where 
Al is the amplitude of primordial gravitational waves, A| 
is the amplitude of curvature perturbations, and we choose 
a pivot point = 0.002 Mpc"'. The tightest published up- 
per limit is r < 0.22 at 95% confidence and is derived from 
a combination of the Wmap five-year measurements of the 
TT power spectrum at low £ with measurements of T ype la 
supernovae and baryon acoustic oscillations (iKomatsu^t al.l 
I2DD9) . 

As a method for constraining r, a direct measurement of 
the BB spectrum has two advantages. First, measurements of 
TT are ultimately limited by cosmic variance at large angu- 
lar scales, and the temperature data from Wmap have already 
reached that limit. Second, r constraints from TT are limited 
by parameter degeneracies; in particular, there is a strong de- 
generacy with the scalar spectral index n,. The inflationary 
BB spectrum, in contrast, suffers little from parameter degen- 
eracies at the angular scales that BiCEP probes, and the BB 
amplitude depends primarily on r. 
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Fig. 1 1. — BiCEP's combined power spectra (black points) are in excellent agreement with a ACDM model (gray lines) derived from Wmap five-year data. 
The (for nine degrees of freedom) and PTE values from a comparison of the data with the model are listed in the plots. The asterisks denote theoretical band 
power expectation values. Power spectrum results from the alternate analysis pipeline are shown by the open circles and are offset in I for clarity. 




Fig. 12. — BiCEP measures EE polarization (black points) with high S/N 
at degree angular scales. The BB spectrum (open circles) is overplotted and 
is consistent with zero. Theoretical ACDM spectra (with r = 0.1) are shown 
for comparison. 



To determine a constraint on r from BiCEP's BB spectrum, 

we compare the measured band powers, '^^j^^, to a template 
BB curve and vary the amplitude of the inflationary compo- 
nent, assuming that it simply scales with r. The template BB 
curves are calculated using fixed ACDM parameters derived 
from Wmap five-year data and a wide range of trial r val- 
ues. We include a constant BB component from gravitational 



lensing, although the contribution is negligible at low multi- 
poles. From the template BB curves, we compute the expected 
band powers, "if **(r). We apply offset-lognormal transforma- 
tions to the data, expected band powers, and covariance ma- 
trix (Equations 25 - 27]i, and we calculate 



-Z*^(r)]^[D^^(r)]- 



-T?\r)-\ 



at each trial r value. The likelihood is then 



£ oc 



1 



V'det[D««(r)] 



-X ('■)/2 



(28) 



(29) 



Figure 14 (left panel) shows that the maximum likelihood 
value obtained from BiCEP data is r = Qffl^^^^. For com- 
parison, the central panel shows the maximum likelihood r 
values from 500 signal-plus-noise simulations with r = input 
spectra. We calculate the upper limit on r by integrating the 
positive portion of the likelihood up to the 95% point, and we 
find that the BiCEP BB spectrum constrains r < 0.72 at 95% 
confidence. This constraint is consistent with those obtained 
from simulations (Figure [14] right panel). 

We have cross-checkeolhe BiCEP r constraint with two 
methods. First, the alternate analysis pipeline yields very sim- 
ilar estimates of r, adding confidence to our result. Second, 
we have generated 100 signal-plus-noise simulations using in- 
put spectra corresponding to r = 0.72, and we have calculated 
r likelihood curves for each of these realizations. Figure 
shows the histogram of maximum likelihood r values, whlcl 
peaks as expected around the input value of 0.72, confirm 
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Fig . 13. — BiCEP's TE, EE, and BB power spectra c omplement existing data from other CMB polarization experiments ("Leitch et al. 2005 Montr oy et al.| 
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BiCEP's constraints on BB are the most powerful to date. 



ing our pipeline recovers an unbiased estimate of r. Only two 
of 100 simulations have maximum likelihood r values that fall 
below the data; this simple alternative statistic to the Bayesian 
95% upper limit suggests at a similar level of confidence that 
the BiCEP data excludes the r = 0.72 hypothesis. 

While limits on r derived from CMB data are still driven 
by the Wmap measurement of the TT power spectrum, Bi- 
CEP's limits on the amplitude of the BB spectrum are an order 
of magnitude more powerful than any previous measurement. 
The improvement in the power of BB to constrain r is illus- 
trated by repeating the above analysis using Wmap BB data. 



where we obtain a limit of r < 6, versus the BlCEP constraint 
of r< 0.72. 

14. CONCLUSIONS 

Motivated by the exciting possibility of detecting, albeit in- 
directly, the gravitational wave background due to inflation, 
many efforts are underway to develop the instrumentation and 
methods necessary to search for the B-mode component of 
CMB polarization at degree angular scales. In this paper, we 
have presented initial results from BiCEP, the first experiment 
optimized specifically to search for the inflationary B-mode 
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Fig. 15. — Maximum likelihood r values are shown for 100 signal-plus- 
noise simulations with input spectra for r = 0.72, a value chosen to match 
BiCEP's 95% upper limit. As expected, the distribution peaks around 0.72, 
showing the BiCEP pipeline recovers an unbiased estimate of r. Only two 
of the 100 simulations fall below r = 0.02, the maximum likelihood value 
derived from BiCEP data, offering an alternative demonstration that r = 0.72 
is strongly disfavored by the BiCEP data. 



combination of these results from BiCEP and of results from 
its sister experiment, QUaD, confirm with significant preci- 
sion the theoretical prediction for the shape and amplitude of 
the ii-mode spectrum on angular scales that span its first six 
peaks. In addition, both experiments have shown that there 
is no B-mode component comparable in amplitude to the E- 
mode component. 

Because BlCEP is designed specifically to probe the de- 
gree angular scales at which the inflationary Z?-mode signature 
peaks, BiCEP's upper limits on the Z?-mode angular power 
spectrum provide the first meaningful constraint on the infla- 
tionary gravitational wave background to come directly from 
CMB B-mode polarization, r = 0.02!jj;^^ or r < 0.72 at 95% 
confidence. Because this limit is dominated by statistical 
noise, constraints on r may be expected to improve in propor- 
tion to observing time (not its square root), and it is likely that 
analysis of the remaining data from the three seasons of ob- 
serving will provide improved limits on r. Future BiCEP anal- 
yses will include the full data set, explore relaxing the conser- 
vative cuts employed for the initial analysis, and optimize the 
power spectrum recovery, all of which will strengthen BiCEP 
constraints on r significantly. 



signal. Analysis of a subset of the first two years of obser- 
vations provides the most sensitive measurement to date of 
CMB polarization at degree angular scales. The TT, TE, and 
EE angular power spectra are measured with high S/N; the 
TB, BB, and EB spectra are consistent with zero. The spectra 
are consistent with a ACDM model using parameters derived 
from Wmap five-year data. The polarization data pass all of 
the statistical tests that we have been able to devise. Further- 
more, the results have been cross-checked by two independent 
analysis pipelines, whose agreement is excellent at all stages 
of the data processing, including maps, power spectra, jack- 
knives, and comparisons with theoretical models. A detailed 
study of BiCEP's potential systematic errors indicate that all 
systematic effects are below the level of the statistical errors 
and are not limiting factors in the analysis. 

BiCEP has detected, for the first time, the first peak in the 
angular power spectrum of the Zi-mode polarization. The 
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